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Abstract 

The paper generalizes the results on soft - hard decomposition of the character- 
istics of QCD jets obtained in by taking into account the effects of fermions and 
running coupling constant. 
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Physics of QCD jets is one of the most actively developing topics in application of 
Quantum Chromodynamics to phenomenology of multiparticle production. The philoso- 
phy and practice of corresponding analysis is described in j^, see also the reviews [SI HUH]- 
One of the central issues in QCD applications to multiparticle production is distinguish- 
ing between the soft nonperturbative and (semi)hard perturbative mechanisms of particle 
production [HI [7|. Soft gluons play a very special role in applications of QCD parton 
model to multiparticle production. This was a main motivation behind the analysis of 
[T], where dynamics corresponding to perturbative evolution of QCD jets was separated 
into soft and hard components according to the energy carried by an offspring parton. 
Such decomposition allows to make a separate analysis of multiplicity distributions arising 
from soft and hard kinematical subspaces and suggest experimental observables allowing 
to verify the theoretical picture. In P only gluon contribution was considered. In ad- 
dition, a simplifying assumption of a fixed (energy independent) QCD coupling constant 
was made. In the present work we lift these two assumptions and study the soft-hard de- 
composition of QCD jets in a full picture that includes fermions and the effects of running 
coupling constant. 

The central object allowing to give a concise and parsimonious description of the 
evolution of multiplicities of QCD jets is a generating function. It is constructed from the 
probabilities of creating n partons in a jet P„: 

oo 

G{z) = Y.p^{i+zr. (1) 

n=0 

Introducing generating function allows to write a transparent system of integro-differential 
equations describing the evolution of QCD jets: 

— ^ = / dxK^{x)^l\GG{y + \nx)GG{y + Hl-x))-GG{y)] 
ay Jo 

+ rif dxK^{x)^l[GF{y + lnx)Gp(i/ + ln(l - x)) - G'g(z/)], (2) 

dG 

— ^ = / dxKf{x)^l[GG{y + \nx)Gp{y + Hl-x))-GF{y)l (3) 
dy Jo 

where y = \n{p6/Qo) is a parameter of jet evolution, p is initial energy, 6' is a jet opening 
angle, Qo = const, as is a strong coupling constant, n/ is a number of quark flavors. 



"^^^^^^ '-^'^ 



^^_ llN,-2Nf _ 17iV,^-iV/(5iV, + 3C^) ^ 

The kernels of the evolution equations (0) and Q read 

K^{x) = - - {1 - x)[2 - x{l - x)], 



1 X 

— 1 + - 

X 2, 



(6) 
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where x and 1 — x are the shares of energy carried by the created partons, Ai'c = 3 is a 
number of colors and Cp = {N^ - l)/2N^ = 4/3. 

The main idea of pQ is to separate soft and hard contributions to jet multiphcity 
{n{y)) = {n^{y)) + {n^{y)) by introducing a cutoff xq <^ 1 in integration over the phase 
space in the evolution equations (0) and Q. The soft contribution (n^) is that coming 
from integration in the interval [0, xq] and the hard one (n^) - in the interval [xq, 1]. Be- 
cause of complete equivalence of two gluons emitted by their parent gluon, the kernels of 
the equations are symmetrical for the substitution of x by 1— x in the initial formulation of 
the problem. It is, however, possible to use this symmetry (see and get the unsymmet- 
rical kernels © at the expense of labeling one of the gluons as an offspring gluon. Then 
the contributions {n^) and {n^) can be easily separated. They differ in such formulation 
and correspond to mean multiplicities of sets of jets with energies less than xop and larger 
than xop. There is no such problem for jets in two-jet events of e+e~-annihilation because 
their energies are fixed from the very beginning. However, subjets in these events, jets in 
three-jet events or jets in hadronic reactions vary in energies. Therefore, it is important 
to know the energy behavior of soft and hard sets of such jets. The above proposal solves 
this problem. 

In [1^ only gluon evolution at fixed coupling constant was considered. Our aim in the 
present paper is to generalize the calculations of P by taking into account the contribu- 
tions of fermions and running coupling constant. The calculations we have to perform 
mirror those in jHl IHI • The only difference is the above-mentioned restriction in the inte- 
gration over the phase space in computing the soft contribution to jet multiplicity. Let us 
remind that a convenient way of taking into account the effects of running coupling fTU| 
in the equations describing the evolution of the moments of jet multiplicity distributions 
with energy is to Taylor expand the generating functions G in the evolution equations in 
perform integration over x and collect all contributions of the same order in the right 
hand- and left hand- sides of the resulting equations. Performing, in complete analogy 
with jHlEl, this calculation for the average soft multiplicities (n^) and {rip) we get the 
following expressions for their second derivatives: 

= lo[{nG)xo + 2{u{nG) - {nG)')wi - {2u {no)' - {nG)")w2 + ^^w^ 
' [{u{nG) - {ug)' - 2u{nF) + 2{nF)')w4, + 2{u{nG)' - 3u{nF)' + {nF)")w^ + {uf)'" Wq]], 



ANc 

= 'll[{nG)xQ + {u {ug) - {nG)')w7 + {2u {ug)' - {nG)")ws - ^{nG)"'wg 

+ (m («g)' + u (uf)' - {nF)")wio - ^{nF)"'wu], (7) 

where u = 2i?7g, B = Pq/SNg, Pq = {llNc — 2Nf)/3, Wi{xo) . . . Wii(xo) are the functions 
of the cutoff Xq listed in the Appendix. 

The multiplicities of soft jets {tiq) and (rip) can be obtained as the perturbative 
expansions in powers of 70: 

(ng) = {ug) [xo - 6i7o - &27o - hlo] , (8) 
= {nF)[xo - rfi7o - (i27o " ^37o]- (9) 
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7 ={nG)'/{nG) =7o(l-ai7o-a27o-a37o) + 0(7o)> 
r = {nc)/ {np) = ro(l - r^o - r27o " ^^370) + ^(70), 



(10) 

fin 



0- 



ro = Cy/C^ = 9/4,7;^-57, 

The numerical values of Oj, are given in [Hj. Using these relations one easily gets the 
second derivative of {tiq) as 



ML 

Mil 



Xo-(&i + (2ai+5)xo)7o-[&2-(35+2ai)6i-(2ai5+a?-2a2)xo]7o+0(7o) (12) 



and the similar expression for {np)" / {np) with bi replaced by di. The terms of the same 
order of 70 must be equal on both sides of the equations. From this requirement one gets 
coefficients bi{xo) and di{xo). 



bi{xo) = 2wi(xo) + -r^— — -W4{xo) - (2ai + B)xo, 

b2{xo) = - (s{ai + BY + 202) xq + 2(ai + B)wi{xq) - W2{xo) + 
Uf 2(ai + 5 + ri) ^5(2:0)1 

(ii(xo) = w^{xQ) - (2ai + B)xo, 

wio(xo) 



AN, 



d2{xo) = -(3(ai + B) + 202)2:0 + (oi + 5 - ri)wr{xo) + Wsixo) + 



(13) 
(14) 

(15) 
(16) 
(17) 



The numerical values of bi, di at Nj = 3 for different xq are shown in Table 1: 

Table 1 



Xo 


bi 


b2 


di 


d2 


0.1 


0.095 


-0.417 


0.004 


-0.42 


0.2 


0.03 


-0.447 


0.003 


-0.748 


0.3 


-0.181 


-0.53 


-0.014 


-1.276 



It is seen that the role of the second correction increases at lower values of xq- The 
expressions for the terms of the third order of 70 are rather lengthy and we do not show 
them here. The results for gluodynamics are easily obtained if Np is put equal to zero in 
expressions for 6j. 

The soft jet multiplicities behave asymptotically as the total multiplicities. Some 
perturbative corrections in powers of y^^^'^ (or 70) appear at lower energies. Soft multi- 
plicities are generally proportional to their softness parameter xq at small Xo with terms 
like Xolnxo appearing in some /j. All these conclusions can be tested in experiment. 

Let us note that in the fixed coupling case considered in P the soft jet multiplicities are 
proportional to the total multiplicity at all energies. The proportionality factor depends 
only on xo- It implies that the dependence of their ratio at lower energies demonstrated 
in Fig. 1 is due to the running property of the coupling strength in QCD. 

This work has been supported in part by the RFBR grants 04-02-16445, 04-02-16880, 
06-02-17051. 
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Appendix 



-^0 



tyi(a^o) = / dxV[x)^ / + x - ^ - -Xq + 

Jo Jo Z Z ^ 

w^ixa) = r dx^^^^^—^-2Vix)hixil-x)l 
Jo X 

wsixo) = / ^ - 2l^(x) (In^ x + In^fl - a;))], 

2 

W5(a;o) = XoQnxo - 1) - a;o(lnxo - ^) + 2a;o(^^ - 

2 1 1 

We{xo) = —xl{9\n'^Xo - Slnxo + 2) - yQn^Xo - -Inxo + -) + XoQn^Xo - 21nxo + 2), 



wvixo) = / dx^{x)^ (1 - = xo - -p, 
io io 2 4 

W8(xo) = XoQnxo - 1) - ^(l^^o - 



ty9(xo) = a;o(ln^xo — 21nxo + 2) — ^[2 In^ — 2 Inxo + 1], 

8 

Wio{xo) = / ]ln(l-a;), 

Jo X 

f^" 1 2 

■J^iifa^o) = / [$(a;) ] In (1 — x) , 

Jo X 

(18) 

The functions V{x) and $(x) are defined in wijWt. 
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Figure 1: Energy dependence of the ratio of soft to total multiplicities in gluon jet 
{ng) / {uq) for various cutoffs Xq = 0.1 (solid line, red) Xq = 0.2 (dashed line, blue) and 
Xq = 0.4 (dotted line, brown). 
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